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cj ! Abstract 

I . We study the mutual interaction between the dark sectors (dark matter and dark energy) of the 

^ i^. Universe by resorting to the extended thermodynamics of irreversible processes and constrain the 

former with supernova type la data. As a byproduct, the present dark matter temperature results 

^ ■ in good agreement with independent estimates of the temperature of the gas of sterile neutrinos. 
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It is widely held by now that the accelerated expansion of the Universe can be traced to 
some invisible agent (possibly the quantum vacuum or some or other scalar field), dubbed 
dark energy (DE), which takes slightly over 70% of the energy budget and is endowed with a 
high negative pressure -see |l| for reviews. The possibility that this agent interacts with the 
remaining fields of the standard model is more natural and general than otherwise whence 
it is receiving growing attention -see [2] and references therein. However, local gravity 
experiments set severe bounds on the interaction with conventional matter (e.g., baryons) 
[3| but there is nothing against a possible interaction with dark matter (DM). Unless there 
exists an underlying symmetry that would set the interaction (coupling) to zero (such a 
symmetry is still to be discovered) there is no a priori reason to toss it away. Moreover, 
the existence of this coupling is undeniable on physical grounds: since DE energy gravitates 
it must be accreted by massive compact objects like black holes and neutron stars. In a 
cosmological context this energy transfer from DE to DM must be small but non- vanishing. 

The coupling between DE and dark matter (DM), first introduced to lower down the huge 



value of t 
problem 



le cosmological constant J4], can provide a mechanism to alleviate the coincidence 
6|, lZI or even solve it 



. Furthermore, it has been argued that an appropriate 
interaction between DE and DM can infiuence the perturbation dynamics and affect the 
lowest multipoles of the CMB spectrum [9|, [lO|. Recently, it has been shown that such a 
coupling could be inferred from the expansion history of the Universe, as manifested in the 
supernova data together with CMB and large-scale structure [uJ, |l2, ll3|, [iJ]. Signatures 
of the interaction between DE and DM in the dynamics of galaxy clusters has also been 
analyzed [ISl . [l6[ . Further discussions on the interaction between dark sectors can be found 

in yy. 

The aforesaid interaction has also been considered from a thermodynamical perspective 
-see e.g. [l9|, |20[. Assuming that the DE can be treated as a fluid with a well defined 
temperature, it has been argued that if at present there exists a transfer of energy between 



DE and DM, it must be such that the latter gains energy_ 



other way around for the second law of thermodynamics 21 1 



rom the former and not the 
to be fulfilled. The authors 



of Ref. [20(1 considered a system composed of two subsystems (DM and DE) at different 



temperatures. In virtue of the extensive property, the entropy of the whole system is the 
sum of the entropies of the individual subsystems which (being equilibrium entropies) are 
just functions of the energies of DE and DM even during the energy transfer process. 

The aim of this Letter is to look more closely into the aforementioned transfer between 
both subsystems by resorting the thermodynamics of irreversible processes as formulated by 



Jou et al. 22]. There, to cope with nonequilibrium (i.e., irreversible) situations, the entropy 
is generalized by allowing nonequilibrium quantities to enter it. Recall that, by definition, 
nonequilibrium quantities vanish at equilibrium and, consequently, this generalized entropy 
reduces to the conventional equilibrium entropy in that limit. 

If DE and DM conserve separately in an expanding Friedmann-Lemaitre-Robertson- 



Walker universe we would have 



Pm + 3if (p„ + p„) = , (1) 

p, + 3i7(l + w;,)p, = 0, (2) 



where the equation of state of DM can be approximately written in parametric form as 23] 



Pm — "^m -/W + 2 '^m J-m ; Pm — '^m J-m y^^B — J-j ("jj 

provided that T^ <^ M. Therefore, pm ~ a'^ and px oc eyi\i j —?>{1 + w^) da/ a. Here, 
Wx = Vxl Px < is the equation of state parameter of dark energy. 
The dependence of both temperatures on the scale factor 

Tm OC a~^ , Tx oc exp / —^Wx da/a (4) 

follows from integrating the evolution equation, T/T = —3H{dp/dp)n for each subsystem. 
The latter is a consequence of Gibbs' equation, TdS = d{p/n) +pd{l/n) and the condition 
for dS to be a differential expression. Since Wx < 0, equations (jl]) suggest that currently 

When the DE and DM components interact with each other, Eqs. ([T]) and ([2]) generalize 



to 

pm + 'iH{pm + Pm) = r , (5) 

and 
p^ + 'iH{l + w^)p^ = -T , (6) 

respectively, where F denotes the interaction term. For F > the energy proceeds from 
dark energy to dark matter. 

For simphcity we assume the specific couphng 

r = 3HXp,, (7) 

with A a small, dimensionless, positive quantity. This kind of models has been considered 



in the literature -see e.g. \14, 124, 



3.g. y 



25|- and show compatibility with observation [26 1. 
Upon neglecting the DM pressure, the ratio between the energy densities, r = Pm/Px, is 
seen to obey 

— = 3{X + Xr + w^r), (8) 

where a; = In a. When Wx and A are constants the latter equation integrates to 

. . _ [r(e,)(A + Wx) + A] e3(^+"'^)(---(..)) - A 
''^'''~ X + Wx ' ^^ 

The subscript (eg) indicates the value taken by the corresponding quantity when DE and 
DM are in thermal equilibrium. 

In the scaling regime, i.e., when the ratio r stays constant, the coincidence of the two 
dark components is understood [g, [tI]. Outside this regime r may be considered piecewise 
constant, especially for z < 20. In particular, A = —f^ about the present time, where 
ro — 3/7 is the current r value. 

When the interaction is taken into account, temperatures of DE and DM evolve as 

Tx = T^eq) e-3('"-+^)(^-^(e.)) , and T^ = T(,g) -^ e-[2+3(A+^.)](x-x(,,)) ^ ^^q^ 

r(eq) 

respectively. It is immediately seen that both temperatures vary more slowly than in the 
noninteracting case (recall that < A < \wx\), see Fig{TJ 
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FIG. 1: Evolution of the temperatures (in Kelvin degrees) of DM and DE in the non- interacting 
case (dotted lines), and in the interacting case (solid lines). In drawing the figure we have taken 
Wx = —1 and A = 0.3. 

Consider an isolated system composed of two subsystems, one being DE and the other 
DM, satisfying 



T 






~dt 



dE„ 



and Tx 



dSx dQx dE. 



dt 



dt 



dV 
dt dt 



:ii) 



respectively, with E„ 
one has 



PmV , Ex = PxV and V = a^. Since the overall system is isolated. 



-Q 



(12) 



dt dt 

where Q = —a^ T is the energy transfer rate. In the comoving volume, this is nothing but 
the energy conservation law: Px + Pm + 3H{pm + Px + Px) = 0. 

In 20|, the entropy of the whole system depends on the energy densities and volume 



only, and in virtue of the extensive property, it is just the sum of the entropies, Sm and Sx, 
of the subsystems. In equilibrium thermodynamics |21[], irreversible fluxes -such as energy 
transfers- play no part and they do not enter the entropy function which is defined for 
equilibrium states only. However, in nonequilibrium extended thermodynamics such fluxes 



enter the entropy function, 5**, which is more general than 5* as it can be defined also outside 
equilibrium J22(]. In this spirit, we postulate that S* = S*{pm, Px,V,Q)- Accordingly, its 
time variation is given by 

where —T{Q) is defined as the derivative of S* with respect to Q and we have used Eqs. 
(TTTj) . For the sake of simplicity, we suppose T{Q) = AQ, with A a semi-positive definite 
constant. 

Assuming r piecewise constant, the energy transfer rate at present time can be roughly 
determined as 



dQ 

dt 



'^^^"'^^\h + H\<Q, (14) 



,„ 87r(l + ro) 

since the Universe entered the accelerated regime only recently. Thus, in the current accel- 
erating stage, the second law in extended irreversible thermodynamics, dS* /dt > 0, implies 

AidQ/dt) > T-' - T^i . (15) 

In view of the evolution of DE and DM temperatures, this inequality can be satisfied for 
any negative value of A. However, for positive values, there must be an upper bound on A 
to guarantee the second law. 

In virtue of the expressions (fill) and (fTOl) . Eq. (IT5l) can be recast around the present time 

as 

- A'H^ (^ + h\> e^-^^t-') e"^-^^ - f r^eq) {r(eg) e2^(-«) e^ + e-°-i^'(=^) e^-i^}"^ , (16) 

where A' = 567T(eq) A/(8007r) and we have set w^ and tq to —1 and 3/7, respectively. 
Approximating the Hubble function by /J ~ Hq — Hq (1 + go) x, with go the current value of 
the deceleration parameter, integration of the above inequality yields 

(17) 
where f{x',r(eq),X(^eq)) = ''^(eg){(^(eg) " f) c^^^^^'c^ + | e"°-^^'(E9) e^'^^}"^. Using X = — ln(l + z) 



we rewrite the inequality as 

where g{z, Z(^eg),rieg)) = '^(eg){[r(eg) " f] (1 + Hei))"^'^ + ^V^ + f (^ + ^(e9))°-^(l + 2)"^-^}~^ 

With the help of the expressions for the luminosity distance and distance modulus, 

[^ dz 

dL = c [1 + z) / — — - , and fi{z) = 5 log [di) + 25 , 

Jo -nl^J 
one obtains the following inequality for the latter 



fi{z) < 5 log 



C{l + Z) / h{z',A',Zeq,r(^eq)) '' dz' 



25, (19) 



'0 

where c is the speed of light and h{z,A',Z(eq),r(eq)) stands for the right-hand-side of Eq. 

m- 

The independent parameters r(eg), Z(_eq), and A' can be roughly appraised by resorting 
the second law and SN la data. Here we shall use the sample of the recent golden SNe la, 
27| . in a similar way as in [28]. Since we have employed a linear approximation 



compiled in 



in the above integration and focused on the present accelerated era, we will use only low 
redshift data, namely z < 0.5. 

Taking the equality sign in flT^ . setting Hq = 72 Km/s/Mpc, and fitting the distance 
modulus to the supernovae data, the best fit values of the parameters are found to be 
A' = 6.42 X 10-^ r^eq) = 1-09 x 10^ and z^^g) = 5.56 x 10^ with x^ = 73.6178102. We 
remark that zi^eq) has only a small infiuence on the total x^ value when latter lies in the 
range 10^ < Z(^eq) < 10^. The value for A' quoted above is the upper bound set by the second 
law. For values exceeding that one the distance modulus curve (depicted in Fig|2]) would 
run below the SNe la data, i.e., the second law, dS*/dt > 0, would be violated in the region 
above the curve. 



As reasoned above (see also 20|]), after the thermal equilibrium between DM and DE 
is lost (i.e., for z < Z(eq)), a transfer of energy between these sectors will arise and T^ 
will increase more slowly than in the absence of interaction and, correspondingly, Tm will 
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with three parameters: 

A'=6.42x10^ 

r =1.09x10^ 

(eq) 

z, =5.56x10' 

(eq) 
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FIG. 2: The graph sets the border above which the second law of thermodynamics would be 
violated. 

decrease more slowly -see Fig. [T] This does not bring the systems to any new equilibrium, 
but it certainly slows down the rate at which they move away from mutual equilibrium. To 
assess quantitatively the relevance of the interaction appearing in the expression for dS*/dt 
we compute the current value of the ratio between the third and first terms in Eq. (TT3|) . i.e., 



R 



\AQ{dQ/dt)\o 



11a 



mO 



{(r(eg) - |e2^(-9) + I e-°-^^(-^) } 



(20) 



\T-^{dQm/dt)\o 20r(e,) 

where a = A' Hq. Inserting the best fitting values for A', r(^eq) and -2(eg), we get Rmo — 1-4. 
Thus, the term coming from the interaction should not be neglected in the expression for 
the entropy variation. 

We next compare the present DM temperature with and without coupling between the 
dark sectors. In the absence of interaction, in virtue of Eq. (j4j), we would have 

2 ^, ^H) 



T 



mo 



T(eg)(l + 2;, 



eq) 



(21) 



^{eq) 



and the thermal equilibrium between DE and DM should have occurred earlier (at larger 
zi^eq)) than in the interacting case. If we consider T ^ z and use the best fitting value of Z(^eq), 



8 



we infer that T^o ~ 10~^ Kelvin. In the noninteracting case, the current DM temperature 
can be even smaller. In the interacting case, the present temperature of DM can be estimated 

as 

T^, = ^^T,,(l + ;2(,,)fi. (22) 

As noted above Z(^eq) cannot be strictly constrained since x^ shifts only by a factor of about 

10"^ when Z(^eq) varies from 10^ to 10^. For (a,r(eg), ^eg) = (1.64 x 10^,0.8 x 10^, 

3.87 X 10^) — what implies x^ = 73.61782324, and we objain T^, ~ 0.596 Kelvin. This 



value is consistent with the estimations of Hansen et al. [29|] based on simulations of the 
number of satellite galaxies and astronomical upper limit of the sterile neutrinos -one of the 
best promising DM candidates- with temperature 

/ q4pV\^/^ 
T^„ = T^„ n^K^ ^ 0.581 Kelvin (23) 



m 



(see Eq. (1) in Ref. l29|), where T^^ = 1.946 Kelvin, fi^ = 0.3, h = 0.72 and the lower 
bound niw = 550 eV 30 1. 

In summary, we have further employed the second law of thermodjTiamics to study the 
coupling between the dark sectors of the Universe. More specifically, by using the nonequi- 
librium entropy, S*, of extended irreversible thermodynamics alongside the gold SNe la data 
of Riess et al. [27] we roughly determined the redshift, Z(^eq), at which DE and DM were 
momentarily in thermal equilibrium as well as the ratio between their energy densities at 
that moment. In addition, we estimated an upper bound on the phenomenological parame- 
ter A'. We would like to emphasize that in the absence of any mutual interaction between 
dark sectors the DM temperature results extremely low. However, due to the interaction it 



can meet the independent astronomical upper bound on warm DM particles 29|. 

Admittedly, our analysis is just preliminary and presents some limitations. In the first 
place, we did not specify microscopically the process by which the two dark sectors couple 
with one another -this lies far beyond the scope of this work. We only proposed a phe- 
nomenological expression for the interaction (cf. Eq. ([7])). Secondly, we expanded T{Q) 
just to first order in Q. A more ambitious treatment should go to higher orders. How- 
ever, this would introduce additional unknown parameters and complicate the numerical 



analysis, something to be concerned with in view of the scarcity and hmitations of current 
observational data. 

This said, we believe our study may well open the way to fuller investigations of the 
DE/DM interaction in the light of nonequilibrium thermodynamics. 
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